Draft version August 7, 2012 

Preprint typeset using I^T^X style cmulatcapj v. 03/07/07 



A FIRST APPLICATION OF THE ALCOCK-PACZYNSKI TEST TO STACKED COSMIC VOIDS 

P. M. Sutter 1,2,3 ' 4 , Guilhem Lavaux 5,6 , Benjamin D. Wandelt 2,3,1,7 , and David H. Weinberg 4,8 

1 Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801 
2 UPMC Univ Paris 06, UMR7095, Institut d'Astrophysique de Paris, F-75014, Paris, France 
3 CNRS, UMR7095, Institut d'Astrophysique de Paris, F-75014, Paris, France 

4 Center for Cosmology and Astro-Particle Physics, Ohio State University, Columbus, OH 43210 

5 Department of Physics & Astronomy, University of Waterloo, Waterloo, ON, N2L 3G1 Canada 

6 Perimeter Institute for Theoretical Physics, Waterloo, ON, N2L 2Y5, Canada 
7 Department of Astronomy, University of Illinois at Urbana-Champaign, Urbana, IL 61801 
8 Department of Astronomy, Ohio State University, Columbus, OH 43210 

Draft version August 7, 2012 

ABSTRACT 

We report on the first application of the Alcock-Paczynski test to stacked voids in spectroscopic 
galaxy redshift surveys. We use voids from the Sutter et al. (2012) void catalog, which was derived 
from the Sloan Digital Sky Survey Data Release 7 main sample and luminous red galaxy catalogs. 
The construction of that void catalog removes potential shape measurement bias by using a modified 
version of the Z0B0V algorithm and by removing voids near survey boundaries and masks. We apply 
the shape-fitting procedure presented in Lavaux & Wandelt (2012) to ten void stacks out to redshift 
z = 0.36. Combining these measurements, we determine the mean cosmologically induced "stretch" 
of voids in three redshift bins, with la errors of 5-15%. The mean stretch is consistent with unity, 
providing no indication of a distortion induced by peculiar velocities. While the statistical errors 
are too large to detect the Alcock-Paczynski effect over our limited redshift range, this proof-of- 
concept analysis defines procedures that can be applied to larger spectroscopic galaxy surveys at 
higher redshifts to constrain dark energy using the expected statistical isotropy of structures that are 
minimally affected by uncertainties in galaxy velocity bias. 

Subject headings: cosmology: observations, cosmology: large-scale structure of universe, cosmology: 
cosmological parameters, methods: data analysis 



1. INTRODUCTION 

Characterizing the nature and history of dark energy is 
perhaps the greatest challenge in the near future of obser- 
vational cosmology. Many elementary probes now strive 
to distinguish a cosmological constant from alternative 
theories of dynamical dark energy or modified gravity. 
Most probes rely o n "standard candles," such as Type 
la supernovae (e.g., Aldering et al.||2002[), or "standard 
rulers," such as radio galaxy diameters ( |l3aly et al.|2009[) 
or baryon acoustic oscillations (BAO) (e.g., |Eisenstein| 
et^[2005l |Reid et al.|[20T0l |Bhattacharya et al.pOllf 
Mehta et al. 2012). Reviews of dark energy probes, cur- 



rent c onstrain t s and fo recasts for future expe ri ments in- 



clude Linder (2003), | Albrecht et all ( 2006 1 , Frieman 
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eTaL] (24)081 ), anT | Weinberg et al.| ^Ol^ T 



Over 30 years ago Alcock & Pacyznski (1979, hereafter 
AP) proposed an elegant alternative approach based on a 
hypothetical population of idealized spheres. Their key 
insight was that since galaxy spatial positions are in- 
ferred from both their angular positions and redshifts, 
these spheres will appear anisotropic if one adopts an 
incorrect spacetime metric. Specifically, because line-of- 
sight distances scale with the inverse Hubble parameter 
H~ 1 (z) and transverse distances scale with the angular 
diameter distance Da(z), their ratio, or stretch, mea- 
sures the value of the product H(z)Da(z). In practice, 
the AP test requires only statistical isotropy of the ob- 
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served structures, so the test can be implemented with 
measures of quasar, galaxy, or Lyo forest clus tering, or 
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features in the redshifted 21 cm spectrum ( e.g., Hui et al 
McDonald fc Miralda-Escude1|T999l |Eriksen et al. 
Nusser||2005HKim fc Crott||2007[ |Blake et al.||2011[ 



2005; Ki 

Reid et al.||2012| ). In this paper we apply the AP test to 
the Sutter et al.|( 2012 ) catalog of voids in the galax y red- 



shift survey s of the Sloan Digital Sky Survey (SDSS; York 
eFaT][2000|. 



lb date, most applications of the AP test have focused 
on the autocorrelation function or pow er spectrum (e.g., 
Ballinger et al.||1996| |Matsubara||2004[ ). Specifically, the 



clearest detections of the AP effect have been found in 
the t wo-point correlati ons of galaxies in the WiggleZ sur- 
vey ( |Blake et al.|201lj) and the B aryon Oscillation Spec- 
troscopic Survey ( |Reid et aL 2012 ). However, a successful 
application of the AP test requires handling the large sys- 
tematic uncertainties caused by peculiar motions, which 
introduce redshift-space anisotropy that must be disen- 
tangled from the AP effect itself. Uncertainties in the 



pec uliar velocity corr ections limit the Blake et al. (2011 1 
and Reid et al. ( 2012 ) studies to large, quasi-linear scales, 
where the statistical uncertainties are relatively large. 

Cosmic voids provide an attractive alte rnative f or ap 
plying the AP test, as first proposed by |Ryden"l (|1995| 
and discussed extensively by |Lavaux fc Wandelt| ( |2012 
Voids are the large, underdense regions that occupy a 
large fraction of the volume of the Universe and are 
a natural consequence of the hierarchical growth of 



structure (Hausman et al. 1983 Thompson & Gregory 
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shapes 1 


Ryden & Melott 1996 


Maeda et al. 12011 Lavaux 


& Wane 


clt 2012), voids avoid the regions oi high velocity 



space correlation function and power spectrum. Indeed, 
modeling of peculiar velocities in voids is particularly 
straightforward since they are still in the quasi-linear 
regime. In addition, the scale of voids is fairly small, with 
typical comoving radii ~ 10 ft, Mpc in a densely sam- 
pled survey, and they have a large filling factor, ampli- 
fying their statistical power relative to other techniques. 
We can therefore measure the mean void shape with high 
precis ion in a large volume survey. |Lavaux fc Wandelt] 



(2012) showed that a statistics-limited void AP test can 
dramatically improve the dark energy constraints from 



the redshift surv ey planned for the Euclid satellite ( Lau 
reijs et al.|2011 1 ; the AP test outperforms the B AO con 
straints from the same survey because the scale of the 
voids is so much smaller than the BAO scale, yielding 
correspondingly more precise measurements. 

We take the void s ample for this ana lysis from the 
catalog described in (Sutter et al. 2012). That work 
constructed a void catalog from the mam galaxy red- 
shift survey ( Strauss et al.||2002| and the luminous red 
galaxy (LRG) redshift survey ( [Elsenstein et al.|2~001[) of 



the SDSS Seventh Data Release (DRY; Abazajian etal 
2009|. We identify voids using a modified versio n ot 



the Voronoi-based ZOBOV algorithm ( |Neyrinck||2008[ ) . To 

compensate for the significant Poisson sampl ing noise 



in shape measurements of individual voids ( Shoji & 
Lee 2012), we instead measure the mean void shape by 



"stacking" the galaxy distributions of our identified voids 
in bins of redshift and radius. The SDSS LRG sur- 
vey is sparse, so at z > 0.2 we can only identify large 
voids, which are limited in number. The combination of 
moderate redshift leverage and limited statistics prevents 
us from making a secure detection of the AP effect in 
this sample, but our proof-of-concept analysis addresses 
many of the practical issues that will also arise in future 
data sets at higher redshift. 

In the following section we give a brief overview of 
our method for measuring void shapes and the applying 
AP te st. We review the properties of the |Sutter et al.| 
(2012) void catalog in Section [3] followed by a presenta- 
tion of the stacked voids in Section |4] We estimate the 
uncertainty in the stretch measurement and present the 
AP test as applied to our void stacks in Section [5] Fi- 
nally, we offer concluding remarks and a brief discussion 
of prospects for future surveys in Section [6j 

2. MEASURING VOID STRETCH & THE AP TEST 

Our defi nitions and procedures clos ely follow those de- 
scribed by |Lavaux fc Wandelt] ( |2012[ ) , who present tests 
on N-b ody simulations and forecasts for surveys such as 
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BOSS (|Dawson et al. 2012) and Euclid (jLaureijs et al. 
Given a galaxy s s 



<y latitude 0, sky longitude 
we transform to a hybrid coordinate 



cz 



= — cos (j) cos 9, 
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sin 6, 



(1) 



where c is the speed of light and Ho is the Hubble param- 
eter at redshift z — 0. Note that our coordinate transfor- 
mation preserves relative distances; in effect, we are per- 
forming a slightly modified version of the AP tes t where 
we measure shapes directly in redshift space (see |Ryden| 
1995 for a discussion). Since the AP test only applies to 



two dimensions (the extent along the line of sight and 
the projected angular extent), we may project positions 
within the void onto a plane: 



d v ■ 



(2) 



Zrel 



where (x rc i, y re i, z IC \) are the galaxy coordinates relative 
to the void barycenter ~X. V : 

x rel = x' - X t ,. (3) 

When stacking voids we place all void barycenters at a 
common point and rotate the galaxies within each void 
about a specified axis so that they all share a common 
line of sight. We then pixelize the density using 10 bins 
within the maximum void size in that stack, which helps 
smooth spurious density fluctuations. 

For each stacked void, we assume a radial profile with 
form 



n(r) 



A +A 3 



(4) 



where r = |x rc i|, R v is the void radius, n is mean galaxy 
number densi ty, and Ap and A% are free parameters. As 



discussed by (Sutter et al. 2012), due to the sparseness 
with which galaxies sample the underlying density distri- 
bution the profiles seen in observations are steeper than 
those in dark matter simulations, so we use a steeper 
form for a fitted curve. Using this radial profile we fit to 
an ellipse given by 

n(d v ,z v ) = min (n + ((d v /a d f + (z v /a z ) 2 ) 2 , s 

(5) 

where no is the density at the center of the stack, 
and a z are the semi-axes along the angular direction and 
redshift direction, respectively, and n max is a maximum 
density value. Our fitting procedure requires an esti- 
mate of the uncertainty on a per-pixel basis. We assume 
that the fluctuations in each pixel are Gaussian with two 
separate standard deviations depending on the location 
inside the stack: 



a(d,z) = h\f IE ^> Xn{d v ,z v )<n n 



otherwise 



(G) 



This gives us a per-pixel sampling uncertainty within the 
void radius and a fixed standard deviation outside the 
void. We keep the latter uncertainty fixed because the 
regions outside the void carry large statistical weight but 
are largely unimportant to the fit. The factor of 1/d 
accounts for the cylindrical averaging of pixels as we form 
the stack. The values of <7o and ct\ are parameters we 
marginalize over in our fits. 

We truncate our stack at R cut = 3R V . We run a Monte 
Carlo Markov Chain to explore the four parameters of 
Eq.([5|), the two standard deviations of Eq.([6), and their 
uncertainties. The likelihood that we must then explore 
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takes the form 

N d N z 

x 2 
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(n(d» 



+ 21og£r(d i ,a i ) J (7) 



where is the number of pixels in the angular direction, 
and N z is the number of pixels in the redshift direction. 
The values di and Zi are the values of the coordinates d v 
and z v at their respective indices. The exploration of this 
likelihood gives us both the measurement of the void el- 
lipticity and the overal l uncertainty associa t ed wit h each 



stack. The analysis of Lavaux & Wandelt 
that the error bars produced from this met 



(|2012] found 
10 d were con- 
sistent with the level of scatter among independent N- 
body simulations, though we will conclude below that 
they underestimate the errors in our data set. 

We translate these ellipticities into a cosmological mea- 
surement by applying the AP test, in which we measure 
the ratio of the length along the line of sight to the angu- 
lar diameter of each stacked void. We will call this ratio 
the void stretch. What follows is a br ief discussion of the 
stretch as a function of redshift; see Lavaux & Wandelt 



(2012) for a more complete derivation. 

We wish to take ratio of a void length along the line of 
sight 5z v to its projected angular extent Sd v . In the sim- 
ple coordinate system of Eq. pi), the projected angular 
extent is related to the angular extent by 5d v = cz89 / H a . 
The angular extent in turn depends on cosmology via the 
angular diameter distance Da{z): 

69 =J^, (8) 



D A { z y 

where 5r v is the comoving radial extent of the void 



In 



a flat universe, the angular diameter distance is equal to 
the comoving line of sight distance D c (z): 

where E(z) = H(z)/Hq. Combining these gives the ex- 
pression 

The comoving line of sight distance 5z v is also related 
to D c (z), and hence Da{z) in a flat universe, via 

Slv 



5z v = 



(11) 



dDj^/dz 

where 5l v is the comoving distance along the line of sight. 
Taking the derivative allows us to identify 



5z v = — E(z)5l v . 



(12) 



In an isotropic universe a stacked void should have the 
same extent in all directions; thus, its angular extent 
should equal its comoving distance along the line of si ght . 
This allo ws us to assume 51 v — Sr v . Combining Eqs. ( 10 1 
and ( 12 ) above leads to our desired ratio: 



Sz v 
Sd v 



H y D A (z)E(z) 



(13) 



We identify the void stretch, denoted by e v (z) for a void 
at redshift z 7 as 

c Sz v 
H Q Sd v ' 



e v (z) 



(14) 



As we stack voids within redshift bins, we assume that 
the stack provides a measurement of the average stretch 
in that bin, (Sz v /5d v ), and we will compare that to the 
average expected stretching in that bin: 



e v (z) 



1 

~Kz 



zt+Az 



e v (z')dz , 



(15) 



where the given bin runs from redshift Zi to Az. 

Throughout we will assume a flat universe with a cos- 
mological constant, which gives a Hubble equation of 



E(z, w , w a ) = (n m (i + zf + n A ) 



1/2 



(16) 



where £l m and Q\ are, respectively, the present-day mat- 
ter and dark energy densities relative to the critical den- 
sity. 

3. VOID CATALOGS 



We take our void catalog from Sutter et al. (2012) 



which is based on volume- limited samples of the JNew 



York Unive rsity Value-Added Galaxy Catalog (Blanton 
et al. 2005). This catalog cr oss-matches galaxies from 
SPSS ( fAbazajian et al. 2009) with other surveys using 



impro ved photometric calibrations (iPadmanabhan et al. 
2008 



We also use the LRG catalog of |Kazin et al. 
. Table m summarizes the volume-limited samples 
used in this work. Additionally, to improve our statistics 
by using as many voids as possible, we merge the four 
samples within z < 0.2 into two samples: diml+dim2 
and brightl+bright2. For the smallest and largest voids 
in each sample these combinations violate our assump- 
tion that the voids are evenly distributed throughout 
each redshift bin. However, we do not use the very 
largest voids in any case, and the smallest voids have 
non-uniform redshift distributions within each sample, 
s o we find that t his do es not strongly affect our results. 

Sutter et al. (2012) produced void catalog s using a 
modified version of the void finder Z0B0V (Neyrinck 



2008), which maps the dens ity using Voronoi tessella- 



van de Weygaert|[2007[ ) and collects these Voronoi 
into zones and voids using a watershed tech- 



tions 
cells 

nique QPlaten et al.|2007| |Aragon-Calvo et~aL|2010| ). To 
remove any potential shape measurement bias due to the 
presence of the mask, we choose the "central" catalog of 
voids, which are selected such that they could not pos- 
sibly intersect any boundary or mask in the survey for 
any given rotation about their barycenters. 

4. VOID STACKS 

We have many constraints for grouping voids 
into redshift and radius bins. We choose redshift 
bins corresponding to the limits of the diml+dim2, 
bright l+bright2, Irgdim, and Irgbright samples. Within 
each redshift bin, we divide the voids into bins of radius 
with the following objectives: 

1. Sufficient numbers - Within each stack we require 
enough voids to sufficiently smooth the projected 
density and increase the signal-to-noise so that we 
can make reliable measurements. While this num- 
ber is not fixed, we have found empirically that we 
require at least ~15 voids per stack for the subse- 
quent stretch measurement to converge reliably. 
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TABLE 1 
Data samples used in this work. 



Sample Name 


Catalog 


-^r,min 




•^m ax 


Number of Galaxies 


Mean Spacing (h 1 Mpc) 


diml 


NYU VAGC 


-18.9 


0.0 


0.05 


63639 


3 


dim2 


NYU VAGC 


-20.4 


0.05 


0.1 


156266 


5 


bright 1 


NYU VAGC 


-21.35 


0.1 


0.15 


113713 


8 


bright2 


NYU VAGC 


-22.05 


0.15 


0.2 


43340 


13 


lrgdim 


LRGs 


-21.2 


0.16 


0.36 


67567 


21 


lrgbright 


LRGs 


-21.8 


0.36 


0.44 


15212 


38 



2. Multiple radial bins - As many independent mea- 
surements within the same redshift bin as possible 
allows us to account for scatter that can develop in 
individual measurements by taking weighted aver- 
ages of these multiple measurements. 

3. Narrowness - Each stack should have a narrow 
radial width; otherwise, the density profile will 
smooth out to the point that our shape fitting rou- 
tine cannot reliably measure the ellipticity. Also, 
smaller voids tend to have more Poisson noise and 
can severely degrade the measurement when com- 
bined with larger voids. 

4. Even distribution - The voids in each stack should 
be evenly distributed in redshift so that we reliably 
measure the mean ellipticity without bias. 

Given these conditions, we select four stacks from the 
diml+dim2 and bright 1+brightl samples and two from 
the lrgdim sample. We discard the lrgbright sample be- 
cause there are not enough voids to construct reliable 
stacks. For each sample, our first stack begins at the 
mean galaxy separation d = 1 ^ 3 , which we take as 



the smallest resolvable void ( Tikhonov fc Karachent- 
12011|). rJeiow 20 /i _1 Mpc our 



sev 



2006 



Platen et al. 



bins have width 4 h~ Mpc; above this we switch to 
8 h Mpc widths to collect sufficient numbers of voids 
(the exception to this rule is the smallest stack of the 
bright l+bright2 sample, where we extend the width to 
include enough voids). We treat the lrgdim sample 
slightly differently due to its much poorer resolution, 
small number of voids, and wide range of void sizes. For 
this sample we reject the smallest voids and construct 
one stack with width 16 h~ 1 Mpc and one with width 
36 ft. _1 Mpc. We do not include the very largest voids in 
each sample because they are difficult to reliably combine 
with smaller voids. However, our bins contain over 95% 
of the voids in each sample, meaning that we are taking 
almost full advantage of the void information available 
in each catalog. 

Within each stack we rescale the voids to the maxi- 
mum void size in that stack; i.e, we multiply all positions 
relative to the void center by R v /Rv,ma,x- This reduces 
the effects of Poisson scatter within the inner wall of the 
stacked void, improving our shape estimation. Where we 
have combined voids from different samples, we normal- 
ize the profile of each void to the mean number density 
of galaxies in the sample before adding it to the stack. 

In Figures [TJ [2] and [3] we show the one-dimensional 
radial profiles for each stack in the diml+dim2, 
bright l+bright2, and lrgdim samples, respectively. We 
construct these profiles by measuring the density within 
thin spherical shells. We also show the radial profiles of 
individual voids. Due to our rescaling of void sizes, the 



profiles do not asymptote to the mean number density of 
galaxies in the sample. 

We immediately note the steep profiles relative to re- 
sults from dark matter-only simulations. This is ex- 
pected due to the effects of Poisson sampling and the 
biasing of galaxies as tracers of density and justifies our 
choice of a quartic radial profile for shape fitting. These 
plots also highlight the necessity of stacking: attempts 
to measure the ellipticity of individual vo ids would be 
nearly impossible. As in the analysis of |Sutter et al.| 
(2012 1, we see that the stacked voids greatly enhance the 
signal-to-noise and generate qualitatively similar profiles 
across many void sizes and redshift ranges. 

We show the two-dimensional stacks used in our shape- 
fitting analysis in Figures |1J [5j and [6| As in Lavaux 
& Wandelt (2012), to improve the signal-to-noise we 
fold the stacked void about the <i-axis and to allevi- 
ate the effects of Poisson sampling we di scretize the 
density. While Lavaux & Wandelt (2012) used fixed 
2 h~ 1 Mpc bins for this discretization, our void sizes are 
much larger, so we choose to scale the bin size with the 
void radius so that there are always 10 bins within the 
maximum void radius. For the lrgdim sample we widen 
the pixel size to dampen the larger Poisson noise. The 
black lines in each plot are the isodensity contours of our 
best-fit ellipse for each stack using the technique out- 
lined above. The contours stop at the calculated value of 
n max , which is the mean density outside the void wall. 
We also list the calculated ratio 5z v /5d v in the plot. 

Even with our discretization there is still significant 
variation in the density, especially along the walls of the 
voids in the compensation region. This does not signifi- 
cantly affect our shape-fitting, which is most sensitive to 
the inner edges of the walls, due to its assumed steep pro- 
file. While all the void stacks have clearly-defined inner 
walls, which is essential for our shape fitting algorithm, 
the compensation regions are not clearly defined. This is 
caused by our rescaling of void size, which clears out the 
inner void regions and extends the outer sections of the 
walls. 

Poisson noise and the small number of detected voids 
make finding a reliable fit in the lrgdim samples difficult; 
we will see that this is reflected in the larger error bars 
compared to the other samples. 

5. VOID STRETCH 

Although the fitting procedure defined by equa- 
tions (|4|-([7| provides an estimate of the uncertainty in 
each measurement, the method makes rather simple as- 
sumptions about the errors in density, namely that they 
are Gaussian distributed (justifying \ 2 likelihood) with 
the errors given by equation |6j . This procedure appears 
to give reasonable statis ticaferrors in the dark m atter 
simulations analyzed by (Lavaux & Wandelt 2012), but 



■5 



z - 0.0-0.1, R = 8-12 h~ l Mpc 



z = 0.0-0.1, R = 12-16 hr l Mpc 



1.2 



1.0 



0.6 



0.4 



0.2 




z = 0.0-0.1, R = 16-20 ft" 1 Mpc 



z = 0.0-0.1, R = 20-28 h- 1 Mpc 




Fig. 1. — Radial profiles at < z < 0.1 /or i/ie dim+dim2 sample. We show profiles of individual voids in grey and the profile of the 
stacked void in black. The legend gives the number of voids in each stack. We rescale each void to the maximum void size within each 
stack. 



here we are analyzing galaxy catalogs, which are sparser 
and to some degree biased tracers of structure. We have 
therefore developed an empirical method of estimating 
errors by creating "incoherent" void stacks: instead of 
aligning voids within a stack to have a common line of 
sight, we assign a random set of Euler angles. After ro- 
tating each void in this fashion, we align the barycenters 
and stack the voids as usual. Since we have removed any 
information about the line of sight by construction, the 
mean stretch of the coherent stack must be unity, and 
the scatter about this value arises from statistical fluc- 
tuations that include both the intrinsic scatter in void 
shapes and the errors in the shape measurements. We 
create an ensemble of 100 such incoherent stacks and take 
the rms dispersion of the stretch values as our estimate 
of the la error. 

We may construct diagrams of the measur ed versus 
expected void stretch via the identity in Eq. (13). Fig- 
ure [7] shows such a diagram where we collect our stretch 
measurements from each stack of each sample and com- 
pare those to the expected mean stretch in that redshift 
bin as a function of redshift, e^(z) (Eq. 
bars shown represent the la scatter in t 



15 1 . The error 
ic 100 incoher- 



ent stacks. For the expected stretch we assume a fiducial 
cosmology of = 0.27, f^A = 0.73, and ho — 0.71, con- 
sistent with the latest WMAP 7-year r esults combined 
with sup ernovae and BAO observations (Komatsu et al. 



2011) 



There is significant scatter in the diml+dim2 sample; 
this is most likely due to the uneven distribution of voids 
within the redshift range. We also see significant scatter 
in the Irgdim sample due to the large amount of noise 
present in the stacked voids. Unfortunately, the statisti- 
cal errors are too large to detect the expected signal of 
the AP effect, an increase in stretch from 1.02 in the low- 
est redshift bin to 1.07 in the highest redshift bin (thick 
grey bars). 

Peculiar velocities have a small but not negli gible ef- 
fect on average void shapes. In N-body tests, (Lavaux 



& Wandelt 2012| find a mean bias of 1.16, with peculiar 



velocities systematically flattening voids along the line of 
sight and reducing the stretch factor. Thin grey lines in 
Figure [7] show predictions that include this suppression, 
which actually di sagree more with our SPS S measure- 
ments. However, |Lavaux fc Wandelt (2012) considered 
dark matter rather than sparse, biased galaxy tracers, 



G 



z = 0.1-0.2, R = 13-20 h,- 1 Mpc 



z = 0.1-0.2, R = 20-28 h' 1 Mpc 




z - 0.1-0.2, R = 28-36 fe" 1 Mpc 



z - 0.1-0.2, R = 36-44 fe" 1 Mpc 




Fig. 2. — Radial profiles at 0.1 < z < 0.2 for the bright 1+ bright 2 sample. See the caption for Figure[T]for a plot description. 



z = 0.16-0.36, R = 54-70 hr l Mpc 



z = 0.16-0.36, R = 70-106 hr l Mpc 




Fig. 3. — Radial profiles at 0.16 < z < 0.36 /or t/ie lrgdim sample. See the caption for Figure^for a plot descriptii 
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Sample: diml+dim2, z = 0.0-0.1, R = 8-12 h 1 Mpc 



Sample: diml+dim2, z = 0.0-0.1, R = 12-16 h 1 Mpc 
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Sample: diml+dim2, z = 0.0-0.1, R = 16-20 h 1 Mpc 



Sample: diml+dim2, z = 0.0-0.1, R = 20-28 7i 1 Mpc 
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Fig. 4. — Stacked voids at 0.0 < z < 0.1 for the diml+dim2 sample. Shown are the stacked voids in our hybrid coordinate system (Eq.pjl. 
Colors indicate n/n, the number density in that bin relative to the mean number density of the sample. The black lines are contours of 
constant density of our fitted ellipse for each stack, and the text box indicates the measured stretch, or ratio of the length along the line 
of sight (redshift direction) to the angular extent of the fitted ellipse. 
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Sample: brightl + bright2, z = 0.1-0.2, R = 13-20 h 1 Mpc 



Sample: brightl + bright2, z = 0.1-0.2, R = 20-28 h 1 Mpc 
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[Stretch™ 




Sample: brightl + bright2, z = 0.1-0.2, R = 28-36 h 1 Mpc 



0.0 0.2 0.4 0.6 0.8 1.0 





Sample: brightl + bright2, z = 0.1-0.2, R = 36-44 h 1 Mpc 
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Fig. 5. — Stacked voids at 0.1 < z < 0.2 for the brightl+bright2 sample. See the caption for Figure|l]for a plot descriptii 



and the voids in their analysis were mostly smaller than 
the ones in our sample (though overlapping in size). Fur- 
ther theoretical work is needed to predict peculiar veloc- 
ity effects in the regime studied here and evaluate the 
tension with our measurements. 

We also performed the same analysis as above above 
but including all available voids, including truncated 
voids near the surveys edges and masks. This also in- 
cludes all voids which, if rotated, would intersect any 
edges. The survey boundaries preferentially select voids 
that lie parallel to them: thus the mask edges will bias 
our results with an excess of voids parallel to the line 
of sight, while the redshift boundaries will bias our re- 
sults with an excess of voids perpendicular to the line of 
sight. We found that we do not recover a strong positive 



AP signal, and instead find measurements that scatter 
around ~ 0.95. Since the surface area of the spherical cap 
which defines the redshift boundary is two to three times 
greater than the surface area of the cone which defines 
the mask edges, we expect our results to be biased be- 
low unity when including all voids. We also increase the 
scatter in individual measurements due to the inclusion 
of many smaller and less well-resolved voids. However, 
with ~ 30% more voids in each stack we do reduce the 
error bars for each measurement by a factor of roughly 

Figure [8] shows the relative likelihood of f^M values in 
a flat universe with a cosmological constant, given our 
stretch measurements of SDSS voids. We calculate this 
likelihood using the weighted average measurements in 
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Sample: Irgdim, z = 0.16-0.36, R = 54-70 ft -1 Mpc 



Sample: Irgdim, z = 0.16-0.36, R = 70-106 h 1 Mpc 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 











IStretch = 0.73 


















































































■ 


















■ 
















































































■ 


V 


> 


\ 

















1.0 




d v /R v , 



Fig. 6. — Stacked voids at 0.16 < z < 0.36 for the Irgdim sample. See the caption for Figurc[4]for a plot descriptii 
1.25p— i 1 i—i ' ' ' 1 1-Or 



R=8-12 lr l Mpc 
R=12-16 /r 1 Mpc 
R=16-20 ft" 1 Mpc 
R=20-28 ft" 1 Mpc 
R=13-20 ft" 1 Mpc 
R=28-36 ft" 1 Mpc 
R=36-44 ft" 1 Mpc 
R=54-70 ft" 1 Mpc 
R=70-106/r 1 Mpc 
Bin Average 
Theory 



0.3 
Redshift i 



0.5 



Fig. 7. — Vbici stretch as a function of redshift. We show the 
measured void stretch (points with error bars) of each stack for 
all samples versus the expected mean stretch in that redshift bin 
(thick grey horizontal bands ) assuming a cosmolog y consistent with 
other recent observations | |Komatsu et al!]|2011^ . The thin grey 
lines indicate the expected stretch including the systematic offset 
of 16% i nduced by peculiar velocit ies that was found in simula- 
tions by |Lavaux fc Wandelt | \2012\ . We assign each radial bin a 
unique color. Also, the radial bins are ordered left-to-right within 
each sample redshift range. Note that we distribute the individ- 
ual points within the redshift range for clarity of plotting only. 
The black points with error bars indicate the weighted mean of 
the measurements in that redshift range. Error bars indicate lcr 
uncertainty. Note that for clarity we have truncated the mean el- 
lipticity line of the Irgdim sample so that it does not overlap the 
brightl +bright2 range. 

each redshift bin (black points in Figure [7]), assuming a 
Gaussian likelihood function. To allow for the effects of 
peculiar velocities, we marginalize over a constant mul- 
tiplicative bias factor with a uniform prior in the range 
[1.0 — 1. 2] (the same range of val ues assumed in the fore- 
casts of Lavaux & Wandelt 2012 ). A positive detection of 
the AP effect would correspond to a rejection of Hm = 0, 




Fig. 8. — One-dimensional likelihood derived from voids. We plot 
the relative likelihood for values of SIm in a flat-A universe, after 
marginalizing over a constant peculiar velocity distortion factor in 
the range 1 — 1.2, assumed to be independent of redshift. 

since a flat, pure-A universe has constant H(z) = Hq and 
Da{z) = cz/Hn (i.e., it really does have the coordinate 
system of Eq.[2|. As expected from Figure [71 our current 
statistical errors are too large compared to the predicted 
stretch signal to detect the AP effect. 

6. CONCLUSIONS 

We have performed the first application of the Alcock- 
Paczynski test to stacked voids to observational data. We 
applied the AP test by measur ing ellipticities of st acked 
voids using the void catalog of |Sutter et al. (2012). The 
stacking procedure greatly reduces the effects of Poisson 
noise and allows us to reliably apply the shape- fitting al- 
gorithm of |Lavaux fc Wandelt ( 2012 ). By grouping voids 
into multiple radial bins we obtain many independent 
measurements, and by dividing the void catalog into red- 
shift bins we obtain measurements across the full range 
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of the SDSS DR7 main sample and most of the LRG 
catalog. However, the limited number of voids and the 
considerable scatter that remains does not allow us to 
positively identify the AP effect over the redshift range 
probed by these data 



The SDSS-III BOSS survey (Dawson et al. 



2012) 



should be a much more powerful basis for void-based 
AP measurements than the DR7 redshift survey ana- 
lyzed here. First, in the range of redshift overlap with 
our Irgdim sample (z = 0.16 — 0.36), the space den- 
sity of BOSS galaxies is a factor ~ 3 higher, which 
enables identification of smaller (and more numerous) 
voids and more accurate measurement of void density 
distributions. Second, BOSS extends this higher sam- 
pling density out to z ~ 0.65, probing a larger comov- 
ing volume (and hence more voids) and reaching red- 
shifts where the predicted AP signal is larger, with a 
stretch factor ey(z = 0.65) = 1.2 for a flat-A universe 
with Q m = 0-27. Future g round-based surveys like Big- 
BOSS (Schlege l et al.|2011| could extend these studies to 
z w 1, while t he space- based emi ssion line redshi ft sur- 



veys o f Euclid ( [Laureijsetal.pOllI and WFIRST ( |Green 
etalp01l| | win" probe much larger comoving volumes at 
1— 2.5 . The Fisher matrix analysis of |Lavaux fc W an- 



delt (2012) implies that a void-AP analysis of the Euclid 
survey should yield significantly tighter dark energy con- 
straints than the BAO analysis of the same survey, with 
a fact or of ten improvemen t in the Dark Energy Task 
Force ( [Albrecht etaL|2006[ ) Figure of Merit. 

Our initial toray into observational application of this 
approach highlights two important directions for future 
investigation. The first is more detailed study of mea- 
surement and parameter fitting procedures and error es- 
timation techniques. Our fitting methods are c losely 
modeled on those of Lavaux & Wandelt (2012), but 
there may be other approaches that make better use 
of the available information, such as using an empir- 
ical radial profile in place of our adopted parametric 
model, changing the radial range of the fit, or down- 
weighting the fluctuations arising from clustered galax- 
ies at the void boundaries. Alternatively, one could 
avoid profile fitting entirely and instead use anisotropy of 
the "void-galaxy cross-correlation function," analogous 
to the cluster-galaxy cross-correlation but centered at 
density minima instead of density maxima. The second 
direction is a more detailed study of peculiar velocity ef- 
fects on mean void shapes, examining its dependence on 
void size and on the spatial and velocity bias of galaxy 



tracers. The likelihood analysis in Figure [8] allows for an 
overall velocity distortion factor and therefore effectively 
uses just the redshift dependence of the signal in Fig- 
ure [7] to constrain cosmology. The goal for future analy- 
ses should be to apply a theoretically computed velocity 
distortion correction to each void sample in each redshift 
bin and marginalize only over the uncertainty in this cor- 
rection, getting an absolute constraint on H(z)Da(z) at 
each redshift. In the cont ext of halo occupati on distribu- 



tion (HOD) models (e.g., Zehavi et al. 2011), we expect 



galaxies to have the same mean velocities as their parent 
halos on average, but the velocity dispersion of galaxies 
could differ from that of the dark matter. This velocity 
dispersion bias ca n itself be constraine d by redshift-space 

so we expect the 



2006) 



galaxy clustering (Tinker et al. 

residual uncertainty in peculiar velocity corrections to 
void shapes to be small, though it may still be the lim- 
iting systematic in void-based AP analysis. 

The statistical errors of this approach are limited only 
by the size and redshift range of spectroscopic galaxy 
surveys, which are expected to grow dramatically in the 
coming years. Cosmic voids are the converse of galaxy 
clusters; primordial density minima expand and deepen 
to form non-linear structures that fill much of the uni- 
verse and are, in a sense, the most "dark energy domi- 
nated" regions of the cosmos. The mean shapes of these 
regions may ultimately provide powerful clues to the na- 
ture of the dark energy that pervades them. 

ACKNOWLEDGMENTS 

PMS and BDW acknowledge support from NSF Grant 
AST-0908902. GL acknowledges support from CITA Na- 
tional Fellowship and financial support from the Gov- 
ernment of Canada Post-Doctoral Research Fellowship. 
Research at Perimeter Institute is supported by the Gov- 
ernment of Canada through Industry Canada and by the 
Province of Ontario through the Ministry of Research 
and Innovation. DW acknowledges support from NSF 
Grant AST- 1009505 and the hospitality of the Institut 
d'Astrophysique de Paris. 

Funding for the Sloan Digital Sky Survey (SDSS) was 
provided by the Alfred P. Sloan Foundation, the Par- 
ticipating Institutions, the National Aeronautics and 
Space Administration, the National Science Foundation, 
the U.S. Department of Energy, the Japanese Monbuk- 
agakusho, and the Max Planck Society. The SDSS Web 
site is http://www.sdss.org/ 



REFERENCES 



Abazajian, K. N. et al. 2009, ApJS, 182, 543 

Albrecht, A. et al. 2006, Report of the Dark Energy Task Force, 

arXiv:astro-ph/0609591 
Alcock, C. & Paczynski, B. 1979, Nature, 281, 358 
Aldering, G. et al. 2002, Future Research Direction and Visions for 

Astronomy. Edited by Dressier, 4835, 146 
Amara, A. & Kitching, T. D. 2011, MNRAS, 413, 1505 
Aragon-Calvo, M. A., Platen, E., van de Weygaert, R., & Szalay, 

A. S. 2010, ApJ, 723, 364 
Ballinger, W. E., Peacock, J. A., & Heavens, A. F. 1996, MNRAS, 

282, 877 

Berlind, A. A. & Weinberg, D. H. 2002, ApJ, 575, 587 
Bhattacharya, S., Heitmann, K., White, M., Lukic, Z., Wagner, C, 

k, Habib, S. 2011, ApJ, 732, 122 
Blake, C. et al. 2011, MNRAS, 418, 1725 
Blanton, M. R. et al. 2005, AJ, 129, 2562 



Daly, R. A., Mory, M. P., O'Dea, C. P., Kharb, P., Baum, S., 
Guerra, E. J., & Djorgovski, S. G. 2009, ApJ, 691, 1058 

Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, 
ApJ, 292, 371 

Dawson, K. S. et al. 2012, The Baryon Oscillation Spectroscopic 

Survey of SDSS-III, arXiv: 1208.0022 
Einasto, J. et al. 2011, A&A, 534, A128 
Eisenstein, D. J. et al. 2001, AJ, 122, 2267 
Eisenstein, D. J. et al. 2005, ApJ, 633, 560 

Eriksen, K. A., Marble, A. R., Impey, C. D., Bai, L., & Petry, C. E. 

2005, Observing Dark Energy, 339 
Frieman, J. A., Turner, M. S., h, Huterer, D. 2008, Annual Review 

of A&A, 46, 385 

Green, J. et al. 2011, Wide-Field InfraRed Survey Telescope 

(WFIRST) Interim Report, arXiv: 1108.1374 
Hausman, M. A., Olson, D. W., & Roth, B. D. 1983, ApJ, 270, 351 
Hu, W. & Haiman, Z. 2003, Phys. Rev. D, 68 



11 



Hui, L., Stebbins, A., & Buries, S. 1999, ApJ, 511, L5 

Jones, M. A. & Fry, J. 1998, ApJ 

Kazin, E. A. et al. 2010, ApJ, 710, 1444 

Kim, Y.-R. & Croft, R. A. C. 2007, MNRAS, 374, 535 

Komatsu, E. et al. 2011, ApJS, 192, 18 

Laurcijs, R. et al. 2011, Euclid Definition Study Report, arXiv: 
1110.3193 

Lavaux, G. & Wandelt, B. D. 2012, ApJ, 754, 109 

Linder, E. V. 2003, Mapping the Dark Energy Equation of State, 

arXiv:astro-ph/031 1403 
Maeda, K.-i., Sakai, N., & Triay, R. 2011, JCAP, 2011, 026 
Matsubara, T. 2004, ApJ, 615, 573 

McDonald, P. & Miralda-Escude, J. 1999, ApJ, 518, 24 

Mehta, K. T., Cuesta, A. J., Xu, X., Eisenstein, D. J., & 

Padmanabhan, N. 2012, cprint arXiv:1202.0092 
Neyrinck, M. C. 2008, MNRAS, 386, 2101 
Nusser, A. 2005, MNRAS, 364, 051027041230005 
Padmanabhan, N. et al. 2008, ApJ, 674, 1217 
Park, D. & Lee, J. 2007, ApJ, 665, 96 

Platen, E., van de Weygaert, R., & Jones, B. J. T. 2007, MNRAS, 
380, 551 

Platen, E., van de Weygaert, R., Jones, B. J. T., Vegter, G., & 

Calvo, M. A. A. 2011, MNRAS, 416, 2494 
Reid, B. A. et al. 2010, MNRAS, 404, 60 



Reid, B. A. et al. 2012, ArXiv c-prints: 1203.6641 

Ryden, B. S. 1995, ApJ, 452, 25 

Ryden, B. S. & Melott, A. L. 1996, ApJ, 470, 160 

Schlegel, D. et al. 2011, The BigBOSS Experiment, 

arXiv:1106.1706 
Shoji, M. & Lee, J. 2012, ArXiv e-prints: 1203.0869 
Strauss, M. A. et al. 2002, AJ, 124, 1810 

Sutter, P. M., Lavaux, G., Wandelt, B. D., & Weinberg, D. H. 
2012, A public void catalog from the SDSS DR7 Galaxy Redshift 
Surveys based on the watershed transform, arXiv:1207.2524 
Thompson, L. A. & Gregory, S. A. 2011, ArXiv c-prints: 1109.1268 
Tikhonov, A. V. & Karachcntsev, I. D. 2006, ApJ, 653, 969 
Tinker, J. L., Weinberg, D. H., & Zheng, Z. 2006, MNRAS, 368, 
85 

van de Weygaert, R. 2007, Voronoi Tessellations and the Cosmic 
Web: Spatial Patterns and Clustering across the Universe, arXiv: 
0707.2877 

Weinberg, D. H., Mortonson, M. J., Eisenstein, D. J., Hirata, C, 

Riess, A. G., & Rozo, E. 2012 
York, D. G. et al. 2000, AJ, 120, 1579 
Zehavi, I. et al. 2011, ApJ, 736, 59 



